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Stability properties and mode signature for equilibria of a model of electron temperature gradient 
(ETG) driven turbulence are investigated by Hamiltonian techniques. After deriving the infinite 
families of Casimir invariants, associated with the noncanonical Poisson bracket of the model, a 
sufficient condition for stability is obtained by means of the Energy- Casimir method. Mode signature 
is then investigated for linear motions about homogeneous equilibria. Depending on the sign of the 
equilibrium "translated" pressure gradient, stable equilibria can either be energy stable, i.e. possess 
definite linearized perturbation energy (Hamiltonian), or spectrally stable with the existence of 
negative energy modes (NEMs). The ETG instability is then shown to arise through a Kreih- 
type bifurcation, due to the merging of a positive and a negative energy mode, corresponding to 
two modified drift waves admitted by the system. The Hamiltonian of the linearized system is 
then explicitly transformed into normal form, which unambiguously defines mode signature. In 
particular, the fast mode turns out to always be a positive energy mode (PEM), whereas the energy 
of the slow mode can have either positive or negative sign. 



I. INTRODUCTION 



An important issue for the stability of equilibria of continuous media concerns the existence of negative energy 
modes (NEMs), spectrally stable modes of oscillation of a medium with negative energy. One reason NEMs are 
important is because equilibria with them, although linearly or spectrally stable, can be destabilized by arbitrarily 
small perturbations. For example, if dissipation is added to the dynamics so as to remove energy from a NEM, then 
it can be proven that the mode becomes spectrally unstable. On the intuitive level, dissipation removes energy from 
the already negative energy of the mode, which makes it more negative and increases the amplitude of mode. In 
non-dissipative systems, NEMs can become unstable with the presence of positive energy modes (PEMs) through 
nonlinear coupling. By this means the system can even develop finite-time singularities while conserving the energy 



of the nonlinear system (see, e.g. 



l| and many original references therein). 



10] on streaming 



In plasma physics, the study of NEMs has a long tradition dating to the early work of Sturrock 
instabilities and Greene and Coppi [2] on magnetohydrodynamic (MHD) type dissipative instabilities in confinement 



systems. NEMs have been studied in many plasma contexts; for example, Vlasov-Maxwell dynamics 3 



drift-kinetic 9] theories, wave-wave interaction in the two-stream instability 



ideal MHD in the presence of equilibrium flows 



14 



11 



12 



8|, Maxwell 



19j, magnetic reconnection 



16j . magnetorotational instability 



17 



13]. 



18| . and magnetosonic 



waves in the solar atmosphere 



3- 



2l| that the Hamiltonian 



Because the concept of NEM is intimately related to that of energy, it was proposed in 
framework is the natural one for defining and investigating this phenomenon, contrary to the usual practice in plasma 
physics. Indeed, once the Hamiltonian structure of the model under consideration is known, an unambiguous definition 
of the energy of the system becomes available: the total energy corresponds to the Hamiltonian of the system and 
the energy of the linear dynamics must come in a natural way from the second variation of this nonlinearly conserved 
quantity. Moreover, the normal form theory for linear Hamiltonian systems, provides a clear and systematic way 
for determining the signature of modes in the neighborhood of an equilibrium of the system. Indeed, for systems 
with discrete degrees of freedom, the Hamiltonian of the linearized system can always be cast, for stable modes, 
into the sum of Hamiltonians of decoupled harmonic oscillators, each of which possesses a characteristic frequency 
and a characteristic signature. Namely this signature, which, for each mode, can be positive or negative depending 
on whether the mode provides a positive or negative contribution to the total energy, provides a systematic way 
to identify PEMs or NEMs of the system. Finally, given the existence of the Hamiltonian structure, energy-based 
methods, akin to 'SW of MHD, can be used to obtain sufficient conditions for stability of equilibria or to indicate 
the presence of NEMs. More precisely, an equilibrium has a NEM, if it is spectrally stable but the second variation 
of its free energy functional, evaluated at that equilibrium, has indefinite sign. 

In this paper, we investigate the presence of NEMs and stability properties of a reduced model for electron tem- 
perature gradient (ETG) instabilities, in the Hamiltonian framework. ETG turbulence has been considered as one of 
the mechanisms that enhances anomalous particle and electron thermal fluxes in tokamaks 



22 



23] . The detection of 



NEMs, is therefore important in order to see what potentially unstable modes might lie dormant in the absence of 
dissipation. When destabilized by dissipation, such modes might enhance the anomalous transport across the confin- 
ing magnetic field. We note, however, that the methods applied for the ETG model, can be applied, in principle, to 
any ideal plasma model. 

The model for ETG turbulence considered here, has been previously investigated in Refs. 



n 

24j,|25j], and in Ref.[2j it 



was described how this model possesses a noncanonical Hamiltonian formulation, with a Poisson bracket that turns 
out to be essentially identical to that for reduced MHD 26( (see also [27]). In the present paper we first provide 
further information about the Hamiltonian structure of the model by deriving explicitly its families of Casimir invari- 
ants. We then determine sufficient conditions for stability of generic equilibria by making use of the Energy- Casimir 



method. The investigation of the presence of NEMs is carried out for the case of homogeneous equilibria, for which 
we derive an explicit condition for the presence of NEMs. This condition warns us that NEMs are present if the value 
of the equilibrium pressure gradient lies in a given interval whose end points depend on the perpendicular wave vector 
and the magnetic field curvature. In particular, the length of this interval shrinks to zero as the perpendicular wave 
number goes to infinity. After obtaining the eigenvalues and eigenvectors of the system, we explicitly construct and 
carry out the transformation that puts the Hamiltonian for the ETG model into normal form. 

The paper is organized as follows. In Sec|H]we review the model and its Hamiltonian formulation, and then derive 
the Casimir invariants. In Sec lIIII we apply the Energy-Casimir method and obtain conditions for energy stability. In 
Sec lIVl after reviewing the theory of mode signature for linear Hamiltonian systems we apply it to the ETG model. 
Explicit conditions for the existence of NEMs, their relationship to energy stability and spectral stability conditions 
coming from the dispersion relation are derived, and the normal form transformation is explicitly obtained. Finally, 
we conclude in Sec El 



II. HAMILTONIAN STRUCTURE OF THE ETG MODEL 



The non-dissipative ETG driven turbulence model of 24j, 125J is given by 



d p 



,V 2 cf> + x] 



>rx 



{Vrx, 4>] 



(1) 
(2) 



where two-dimensional slab geometry is assumed so that the field variables, <f> the stream function and p the pressure, 
are functions of the Cartesian coordinates (x,y). The quantity [/, g] = df /dxdg/dy — df /dydg/dx is the canonical 



Poisson bracket. Equations (JTJ) and © are written in the normalized form described in 
parameter r is defined as 

Ll 



241, where the constant 



LbLp ' 



(3) 



with L n , Lb and Lp being the characteristic length scales of variation of the background density, magnetic, and 
electron pressure fields, respectively. The parameter r is thus related to the mechanism providing the drive for the 
ETG instability. In particular, in the limit r — ► 0, which is attained if one flattens the electron pressure gradient, the 



linear dispersion relation indicates that ETG modes degenerate into marginally stable drift waves 



25] 



In 24 1, the authors showed that the system CO)-© possesses a Hamiltonian structure in terms of a noncanonical 



Poisson bracket. This means (see, e.g. 1 



27]) that the system can be cast in the form 



dt 



= {x\H), 



i = 1, • • • , n 



(4) 



with x J ( x ! t) indicating a suitable set of n field variables (with n = 2 in our case) and H [x 1 , • • • , X™] a Hamiltonian 
functional that is conserved by the dynamics. The Poisson bracket {, } appearing in Q is an antisymmetric bilinear 



binary operator satisfying the Leibniz rule and the Jacobi identity. For the model of ([1} and ([2]), it was shown 
that, with the choice x 1 = <P, X 2 — the Hamiltonian of the system is 



24] 



H[4,,p] = ljd 2 x (> + |V0| 2 -^ 



(5) 



and the Poisson bracket is 



{F, G} = - J d 2 x(4> - V 2 - aOI/T 1 ^, C^Gd + (p + rx) ([C^F^ G p ] + [F p , C~^\), 



(6) 



where the operator C, and its inverse £ _1 are formally defined so that Cf = f — V 2 /, and £ _1 £/ = ££ _1 / = /, for 
a generic field /. The subscripts on F and G in ([6]) denote functional derivatives with respect to the fields 4> or p. 

Noncanonical Poisson brackets such as ([6]) are characterized by the presence of so called Casimir invariants (see, 
e.g. Ref. [l(), due to degeneracy in the cosymplectic operator of the bracket. More precisely, a Casimir invariant of a 
Poisson bracket is a functional C(x 1 , ■ ■ • >X n ) that satisfies 



{C,F} = 



(7) 



for any functional F of the field variables. Because they commute in particular with any H, Casimir functionals 
are preserved during the dynamics. In order to derive the Casimir invariants of the ETG model, it is convenient to 
introduce the variables 



V=-?=+Vrx, \ = <j>-V 2 <p~ x, 



(8) 



which correspond to a "translated" pressure and to a variable analogous to the potential vorticity of the Charney- 
Hasegawa-Mima equation 28|, l29( , respectively. In terms of these variables the model equations read 

d 



dt 

d_ 

m 



A = -[C-\X + x),X] + [V,Vrx], 
V = [PX^iX + x)}, 



(9) 
(10) 



whereas the Hamiltonian and the bracket become 



H{W) = \ I d 2 x((\ + x)£- 1 (\ + x)-V 2 + 2^Vx) , (11) 



{F, G} = - J d 2 x\[F x ,G x ] + V ( [Fx , Op] + [F V ,G X ]) , 



(12) 



where the bracket is seen to be identical to that for reduced MHD as first given in 26]. 

Applying (J7J we deduce that the equations determining the Casimir invariants for our system are 

[C x ,X\ + [C v ,V] = 0, 

[Cx,V}=0. 

By solving (fT3)) we see that the system admits two independent infinite families of Casimirs: 



(13) 



Ci = J d 2 xH(V), C 2 = J d 2 xXF(V), (14) 

with H and T arbitrary functions. The dynamics described by the inviscid ETG model is then subject to an infinite 
number of constraints imposed by the conservation of the Casimir invariants (1141) . For instance, as a consequence 
of the conservation of C%, integrals of the potential vorticity A over regions bounded by contour lines of V will be 
conserved during the dynamics (see [301). 

n 

Notice that the constant of motion X found in [24| is given by 



X = J d 2 x(y + 2{(j) - V 2 <t>)p 

= 2y/F J d 2 x\V + J d 2 xV 2 - 2r J d 2 xx - V 2 </> - |) , (15) 



which is a linear combination of two particular Casimirs of (|14l) with the realization that the time derivative of 

-2r J d 2 xx - V 2 - |) (16) 

vanishes if d<j>/dx and dtfi/dy vanish or are periodic at the boundaries. Therefore, X amounts to a single special case 
of the general families of (fl4"|) . 

III. ENERGY STABILITY 



The Hamiltonian formalism provides a systematic procedure for implementing the Energy-Casimir method for 
investigating stability of equilibria (see, e.g. [l, 31, 32 1), a stability method that originated in plasma physics in 33 1 
that has often been referred to as nonlinear stability. This method has been adopted in many works; for example, in 



the context of fluid models for plasmas in [34^37]. It provides sufficient conditions for stability by taking the second 
variation of the free energy, the Hamiltonian plus Casimir invariants, and extracting conditions that are necessary 
for definiteness. Because the method is based on nonlinear constants of motion, the stability conditions obtained are 
stronger than conditions that emerge from dispersion relations, i.e. spectral stability conditions, that follow entirely 
from the linear equations of motion. Indeed this second variation stability, which we will refer to simply as energy 
stability, implies linear stability 



i|, but the converse is not true. In some works energy stability is called formal 
stability when an additional convexity estimate is not provided. Usually these estimates are rather trivial and even 
when they are provided they are only a small part of a mathematically rigorous stability proof - for this reason we 
eschew this terminology. 

For noncanonical Hamiltonian systems, equilibrium solutions can be found by solving the equations that result 
from extremizing the free energy functional. For our system the free energy functional is given by F = H + C\ + C2 
(which is not to be confused with the generic functional F of our Poisson brackets), with H given by (|12[) and C1.2 
by p4[) . For convenience wc introduce the new variables according to the transformation 

K = \ + x, V = V (17) 

and then drop the bar on V in the following. In terms of these variables the Hamiltonian and the bracket become 

H(A,V) = ~ J d 2 x (AC- 1 A - V 2 + 2^/rVx) , 
{F,G} = J d 2 x(x-A)[F A ,G A ] -T([F a ,Gv} + [F v ,Ga}). 

The free energy functional is then explicitly given by 

F(A, 'P) = \J <?x (A£ -1 A - V 2 + 2y/?Vx) + J d 2 xU{V) + J d 2 x(A - x)F(T) , (18) 

and the equilibrium equations, obtained from setting the first variation SF equal to zero, are 

F A = C~ X A + F{V) = 0, (19) 
F v = -V + ^ + n'{V) + {A-x)F'{V)^0 7 (20) 

where the prime denotes derivative with respect to the argument of the function. Due to the presence of the arbitrary 
functions in the Casimirs, such equilibrium equations possess free functions. Specifying these selects from a class of 
equilibrium solutions. In particular, choosing T corresponds to fixing the relation between the equilibrium stream 
function (f> eq = C~ x A eq and the equilibrium translated pressure V eq . 



As indicated above, an equilibrium solution of (fT9|) - ([20|) is energy stable (and therefore linearly) stable, if the 
second variation of F, evaluated at that equilibrium, has a definite sign. In terms of the variables <f> and V , the second 
variation of F reads 



S 2 F = J d 2 x\ i (l-F / (V))\£- 1 5A\ 2 + (l-2F'(T))\£- 1 VSA\ 2 +F'(V)(SA + ST) 
-F'iV^C-^SA) 2 + ( H"(V) - 1 + (A - x)F"(V) - F'(V)) \5P\ 2 



(21) 



From (|2Tjl one immediately obtains sufficient conditions for positive definiteness of 6 F in the case of no flow: 



F(V eq ) = and U"(V eq ) > 1 . 



(22) 



From (fl"9f . F(V eq ) = implies no flow, while W'('Peg) > 1 gives a condition on the pressure profile. Generally 
speaking, the situation with flow, when F(V eq ) ^ 0, is expected to have NEMs 1[. However, this case is more 
complicated to analyze, with the Poincare inequality often being of use, but we will not pursue it further here. 

IV. NEGATIVE ENERGY MODES 



In this section, we first review the theory of NEMs in the finite degree-of-freedom Hamiltonian context 39|, since 
it applies directly to finite systems with discrete spectra (e.g. 



11 



complicated 5 



40-42 



The situation for continuous spectra is more 



and will not be considered here. Subsequently, after carrying out a spectral stability analysis 
of the system linearized around homogeneous equilibria with no flow, we make use of the Hamiltonian formalism in 
order to detect the presence of NEMs among the stable modes of the linearized system. Finally we carry out the 
explicit transformation that casts the corresponding Hamiltonian into normal form. 

A. Review of mode signature and normal forms for linear Hamiltonian systems 



A real canonical Hamiltonian linear system with N degrees of freedom is generated by the canonical Poisson bracket 

Of dg df dg 



if, 9} 



f dpi dpi dq 1 



(23) 



and a quadratic Hamiltonian 



(24) 



8 

where z — (qi, qN,p\, ...,pn), Aij are the elements of a 2N x 2N matrix with constant coefficients, and repeated 
sum notation is assumed with i, j = 1, 2, . . . , N. The resulting equations of motion can then be compactly written as 

z = J c Az, (25) 

where 

is the 2N x 27V canonical symplectic matrix (cosymplectic form). Assuming 

z = Se iuJt + z*e~ iujt , (27) 

with * indicating complex conjugate, (|25[) yields the eigenvalue problem 

iuj a z a = J c Az a , a = 1, • • • ,N, (28) 

where we have dropped the tilde on the eigenvectors and have added an eigenvalue label a. In what follows we 
will assume distinct eigenvalues, precluding the existence of nontrivial Jordan form and possible secular growth in 
time. We also assume that the eigenvalues ui a are real, which is the case of interest for detecting mode signature. 
Because our dynamical variables are real, the remaining N eigenvalues are given by w_ Q = — uj a and the corresponding 
eigenvectors are z a — z* a . Defining fl := J^ 1 , the symplectic two- form, we construct the quantity 

h(a,/3) :— iuj a zjQz a — zjAz a , (29) 

where T denotes transpose. 

It can be easily shown that the property h(a, (3) — h(/3, a) = holds. Then, from this relation and the antisymmetry 
of fl, it follows that 

h(a,j3)=0, if -a. (30) 

On the other hand, 

h{— a, a) — z 1 i a Az a = z^ T Az a — iuj a z^ T Qz a , (31) 

is clearly the energy (Hamiltonian Hl) of the mode (z a , u: a ; z* , — uj a ). Evidently, z* T f2z Q is a purely imaginary 
number, and a normalization constant for the eigenvectors can be chosen in such a way that 

zfriz a = ±2i, (32) 



with the sign, an invariant, depending on the specific mode under consideration. Note that the left-hand side of (|32j) 
represents the Lagrange bracket (symplectic two- form) of z* and z a . If z a is an eigenvector, associated with a positive 
eigenvalue u> a , and 

z*Jnz a = -2i, (33) 

then (z Q , u a ; z* , — uj a ) corresponds to a positive energy mode, otherwise it is a negative energy mode. This can be 
easily seen by observing that, in the case of a PEM, the corresponding energy is given by 

h(— a, a) = iui a z^ T Qz a = 2u a > 0. (34) 

Note that, although here we carried out an analysis with canonical coordinates, Sylvester's theorem guarantees that 
the signature of a mode (i.e., whether it is a PEM or a NEM), does not depend on the choice of the coordinate system. 
The distinction between positive and negative energy modes becomes even more transparent when we are re- 



minded that, for stable modes there exists 3, 



11 



39j a canonical transformation T : (Qi, ■ • • , Qn, Pi, • ■ • , Pn) 



{Qi, • • • In, Pi, • • ~Pn), that casts the quadratic Hamiltonian of (|24p into the following normal form: 

1 N 

H L = -Y,v a u a {P 2 a + Qi) 1 (35) 

a=l 

where uj a represents the positive eigenvalues of the linearized system, whereas <Xj £ { — 1,1} is the signature of the 
mode. 

If the system contains unstable modes, then they have a different normal form. However, if the Hamiltonian is 
restricted to the stable modes, then it can be written as the Hamiltonian for a system of N harmonic oscillators with 
different frequencies. The modes for which a = — 1, which give a negative contribution to the total energy, correspond 
to the NEMs, while those corresponding to cr, = 1 are, of course, PEMs. 

Once the eigenvalues and eigenvectors of the system are known, the procedure for constructing the map T is 
algorithmic. First one needs to select, among the 2N eigenvectors of the system, the N eigenvectors z a that satisfy 

z* a T VLz a = -2i, a = l,---,N. (36) 

Then the 2N x 2N matrix that defines the transformation is given by 

T = col (Re zi, Re z-i ... Re zjv, Im Z\, Imz2, . . . , Imzjv) , (37) 

which is the matrix with columns given by Re Z\ etc. It can be shown that the transformation constructed in this 
way is canonical and indeed provides the desired diagonalization. 



10 



B. Mode signature and stability for the ETG model 



Now consider a special case of the no-flow equilibria of Sec. IIII) the homogeneous equilibria: 

A eq = £- x A e „ = 0, V eq = a v x , (38) 



where a-p is a constant. The equilibrium solution (|38p corresponds to the choices 



2 \ a-p 

for the Casimir functions that appear in equations (fT9|) and ([20)). 

Linearizing the model equations around this equilibrium gives the system 

dy dy 
dy 

Expanding the perturbations as Fourier series, as follows: 



A = A= A k (t)e- 4kx , V = a v x + P = a v x + ^ P k (*K ikx , (40) 



k— — oo k— — oo 

yields the amplitude equations 



A k = i— ^2-A k + iVtk y T k , (41) 

A = -lap-^^Ak, (42) 

whence, the dispersion relation for modes of the form e i( ^ _k ' x ) is obtained, 

^-T^V + »,^f = 0. (43) 



This expression is in agreement with that obtained in Ref. 
and are given by 

k 



251 ] . The eigenvalues correspond to a slow and a fast mode, 



U! 



(l-Jl-4(l + kl)a v V^), (44) 



2(1 + it* ) 

where we have set k y — k. The corresponding eigenvectors are 

1 + k 2 ~ u 1 



2(l + fci)' 

; ~ 9a + fc 2 ) (1 + \Z 1 -^ 1 + k l) a ^), (45) 



l '>t = -<1^t f ^ = -2a^ {1 ± V 1 - ^ + ki)avVt)Kp (46) 



11 



for k > 0. The system also possesses the eigenvalues u)_ a _j = —w^f, whose eigenvectors are the complex conjugates 
of those of (Hp) . 

From (j4"4")) - (l4"5j) we obtain a necessary and sufficient condition for spectral stability viz. 



1 - 4(1 + ki)a P <Jr > 



(47) 



4(l + fc2)VF 

Therefore, if the electron pressure or magnetic field gradients are such that ap < 0, such equilibria are always spectrally 
stable V k. If qj> > 0, on the other hand, the equilibrium will be stable only for k such that < a-p < 1/4(1 + k^y/r 
is satisfied. In other words, there will always be instability for sufficiently large k±. 

The linearized system of (|41[) and (1421) inherits a Hamiltonian formulation from the nonlinear system, one that can 
be written in the framework described in Sec. IIV Al We can then take advantage of this fact in order to see whether 
NEMs are present in the system and to cast the Hamiltonian into its normal form. 
By using the relation (see, e.g. 



11 



43 



SF 

sK 



k—-\-oo 

£ 

k— — oo 



SF 



,-iky 



- y 



dF 



-iky 



(48) 



where F(A) = F(A&), it can be shown that the Hamiltonian structure of (|41[) and (H2J), is given by the bracket 



{F,G} = J2 



k=l 



ik 
2^ 



dF dG 
d~A k dk- k 



dF dG 



f dF dG dF dG dF dG dF dG 
\dA k dV-k + dVkdA- k dV-k dA k dA- k dV k 



(49) 



and the linear Hamiltonian, which is proportional to 8 2 F 



+oo 



(50) 



k=l 



k=l 



where we have suppressed the sum on k x (note that k x only appears in the combination k\ = fc^ + k 2 ). Although not 
canonical, this formulation, in principle is sufficient in order to detect the presence of NEMs for the equilibrium under 
consideration. Indeed, as already done for the four-field model of 



13j , we can make use of the property that NEMs 
(1 

can change their signature only if they become unstable through a "Krem bifurcation" [1[, or if the corresponding 
eigenvalues go through zero frequency. It is then sufficient to identify a NEM in a particular limit, and we are 
then guaranteed that its signature does not change as long as one of the two above mentioned phenomena does not 
occur. Sylvester's theorem also guarantees that the signature is independent on the choice of the coordinates we 
make. If we fix a wave vector k, then we can first evaluate the energy associated with the corresponding mode in 
the (Afc, A_fc, Vk, V-k) coordinates by inserting eigenvalues and eigenvectors associated to k in the expression for H*. 
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0.5 1.0 1.5 2.0 2.5 3.0 



-0.5 
-1.0 

FIG. 1: Depiction of two possible mode signature situations for ETG modes, depending on the sign of a-p. The plot on the left 
refers to the case of negative av {av ~ —0.3, *Jr = 0.2). In this case both modes are PEMs and the system is energy stable. 
In the plot on the right, av is positive («p = 0.5, yfr — 0.2). For < k± < 1.22, the equilibrium is stable, but the slow mode 
is a NEM. A Krem bifurcation occurs at k± — 1.22. 



This results in 



H 



2tt 



1 - A^Fav - 2k\av\fr ± y/l - 4(1 + k\)a Vy /r 



2a 2 v 



l^/l 2 , 



(51) 



which is the energy of the slow and the fast modes of wave number k (summed over k x ). In order to identify PEMs 
and NEMs, it is sufficient to consider the limit k± — > 0, which yields 



tt fe,fei=o ( 1 - 4:y/¥a v ± y/1 - 4ayy/r \ ~ k 2 
H *>.,f = 2 - ^ j P.jl ■ 



(52) 



We can then sec that, for the fast mode (corresponding to the + sign in (|52|) ), HL^ k±=0 is positive, and therefore 
a PEM. Also, it will remain a PEM as parameters are varied in a continuous way, until the instability threshold is 
reached. For the slow mode, two cases exist for finite a-p. If a v < 0, then H L k s ' k± =° > and, again, we have a PEM. 
If < ap < l/(4y/r), on the other hand, the slow mode is a NEM. 

The above mentioned instability, occurring at large fcj_, for ap > 0, indicates a Krein bifurcation, which is one 
of the possible types of bifurcations occurring in Hamiltonian systems. These situations are illustrated in Fig. [T] 
When ujg is real and negative, that is for ap < 0, both the slow and the fast modes are PEMs and the two branches 
correspond to two dispersive waves with opposite sign. These waves correspond to inviscid drift waves modified by 
the presence of the ETG and magnetic field curvature. In the limit r — > 0, corresponding to vanishing ETG, the two 
modes degenerate into a Hasegawa-Mima drift wave. When oj k is real and positive, i.e. for < ap < 1/4(1 + kjjy/r, 
both modes are still stable but the slow mode is now a NEM. Comparing the two plots of FigHJ one observes that 
uj k went from negative to positive, i.e. it crossed through zero frequency, while changing from a PEM to a NEM. 
For the parameters chosen for the figure, the instability threshold, due to the presence of ETG, occurs at k± « 1.22. 
For perpendicular wave numbers above this value, the equilibrium is indeed unstable. The transition of the two real 
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eigenvalues into a complex conjugate pair, occurring at k± 1.22, is an example of a Krem bifurcation. Note that, as 
predicted by Krein's theorem (see, e.g., Ref. l|), if a Krem bifurcation between two eigenvalues occurs, one of them 
must be a NEM. This is indeed our case. 

If we consider now the stability condition (|22|) for no-flow equilibria, together with (l39l) . we find that the homoge- 
neous equilibrium (|38|) is energy stable if H"(V e q) > 1, which is equivalent to ap < 0. Indeed, if this condition is 
satisfied, the equilibrium is stable, with no NEMs. If we now use the actual pressure p, as a variable, and consider the 
corresponding homogeneous equilibrium p eq — a p x — ^/r(ap — \/r)x, we can reformulate our results in the following 
way. If the equilibrium pressure gradient is such that a p < —r (i.e., ap < 0), then the system is energy stable for 
every k. If, on the other hand, —r<a p < — r + 1/4(1 + fcj_), (i.e. < ap < 1/(4(1 + fc^)y / r') then the system becomes 
unstable through a Krein bifurcation at a critical k±, for fixed r. For k± below this critical value, the equilibrium is 
spectrally stable but not energy stable. Indeed, the slow mode in this case is a NEM, which makes the equilibrium 
fragile with respect to the addition of dissipation or nonlinearities. 

C. Normal form for the Hamiltonian of the ETG model 

The explicit knowledge of the eigenvalues of the linearized system, makes it possible to cast the Hamiltonian into 
its normal form, as discussed in Sec. IIV K\ First of all we point out that the transformation (Afc, A_/-, Vk, P-k) ~> 
ill » 9fe , Pi, Pk ) ' corresponding to 



q\ = Jj^{Tk + apA k +V-k + apA- k ), 



Pk = - *i/j£~~2"(^ fe + a vAk. - V-k - apA- k ) 



4 = y^(A fc + A_ fc ), p£ = ty-(A fc -A_ fc ), (53) 

puts the system into canonical Hamiltonian form. Indeed, in terms of the variables {q\, QkiPkiPk)' the bracket of (|49[) 
takes the canonical form 

u~L ^Qi dpi dp\ dqi dq\ dp\ dp\ dq\ ' 
The Hamiltonian (1501) in terms of the new variables, on the other hand, reads 

oo 4 

^=2^^ A ^^' (55) 

k=l 
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where 



A k 



( a c 
c b 
a -c 

V -c 6 



z fe = (gi,«2>Pi>P2) 



(56) 



with a = —^/ra-pk, b = fc(l/(l + fc^) — -y/ra-p), and c = \fr\a-p\k. We emphasize that the transformation (|53[) is 
constructed in such a way that the canonical variables q k and p k are real. We are thus in the framework depicted in 
Sec. IIV Al For each k, the equations of motion are given by 



z k = J c A k z 



(57) 



Upon writing the variables as 



Ql,2 



^l,2 e "t" *1,2 



6-\ O C « 



Pi, 2 — z 3,4 e + z 3,i e ) 



(58) 



and dropping the tilde in the following, we obtain that the eigenvectors of the linearized system are 



7 k - n k 
z s — His 




7 k - n k 
z -s ~ His 



( 1 

i 



7 k - n k 

z f - its 




h h * 

z -f = lif 



( 


1 




—B 




i 


\ 


iB^ 



(59) 



These are the eigenvectors corresponding to the eigenvalues u> k , w_ s , u>f and respectively. In (|59|) we introduced 



the quantities 



B± 



b + a± ^(b + a) 2 - 4c 2 



2c 



(60) 



whereas gi k ^ are complex coefficients. Following Sec. IIV Al the Lagrange bracket for the slow mode reads 



*Vn«; = 2i(l-Bi)gf, q\ 



(61) 



Direct calculation shows that 



V(& + g - 4c 2 ( ^(fr + a) 2 - 4c 2 -(b + a)) 



(62) 



To obtain the inequality in (|62[) . we made use of the fact that, for stable modes, y(&~+ a) 2 — 4c 2 < 6 + a. 

The inequality (|62j) tells us that, for slow modes, we must pick up the plus sign in the general expression ([32 

Moreover, if we choose the normalization constants so that 
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we obtain that the Lagrange bracket for slow modes becomes 



z k _ s T n z k s = 2<. (64) 

Of course, 

zfilzt. = -2*, (65) 

and consequently, according to the definition, given in Sec. IIV A| we have a PEM, when u k _ s = -w k > 0, and a NEM 
when u>_ s = —oj k < 0. This confirms the results we obtained in Sec. IIV 51 with noncanonical variables. 
Following the same procedure for the fast mode, we find 

z k f T Qz k = 2i(l - B 2 + )q k f *q k f . (66) 

Given that 1 — B? < 0, the Lagrange bracket for the fast mode becomes 

z h L f T ttz k f = -2i, (67) 

after having chosen the following normalization for the eigenvectors: 



Because uj k is always positive, according to the definition, (|67p tells us that the fast mode, as expected, is always a 
PEM. 

The transformation that casts the Hamiltonian (|55[) into normal form, following Sec. IIV Al will be a real canonical 
transformation T k that, for each k, maps a new set of coordinates z k — (Q k ,Q2,P k ,P k ), in terms of which the 
Hamiltonian is diagonal, into z k . 

After noticing that 

z k _* T nz k _ s = -2i, (69) 
zf T ttz k f = -2i, (70) 

the matrix associated with the application T k is constructed in the following way: 
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Direct calculation shows that 





\ 


o 



(72) 




\ 




Consequently, the Hamiltonian (|55[) can be finally written as 



+oo 4 



2E E {T k z\ T A%{T k z k ) 3 



k—1 — l 



^ +00 



) w ; (gf + f* 



) 



(73) 



k=l 



The expression (|75|) corresponds to the normal form for the Hamiltonian of the linearized ETG model for stable modes. 
It clearly shows how the corresponding energy can be decomposed into the sum of energies of harmonic oscillators 
which possess, as characteristic frequencies, those of the fast and slow modes. The harmonic oscillators associated 
with the fast modes always provide a positive contribution to the total energy. Those associated to the slow modes, 
on the other hand, give a negative contribution if ui k > 0, which translates into the condition on the equilibrium 
pressure gradient discussed in Sec. IIVBI 



By making use of the Hamiltonian formalism, we have analyzed the mode signature and the stability properties of an 
ETG fluid model. The families of Casimir invariants of the model were obtained, thereby showing that the dynamics 
of the model is subject to an infinite number of constraints. A stability condition has been derived, according to which, 
the absence of equilibrium flow and a restriction on the pressure equilibrium profile imply stability. Subsequently, after 
reviewing the concept of mode signature in the Hamiltonian framework, we have explicitly determined the energies 
of stable modes. From the stability viewpoint, the dispersion relation gives us a spectral stability condition which, 
however, does not give us information about the stronger condition of energy stability. Indeed, our analysis shows that 
spectrally stable homogeneous equilibria can be of two types, depending on the value of the parameters. If a-p < 0, 
equilibria are spectrally and energy stable (i.e. with no NEMs). If < ap < 1/(4(1 + k^_)^/r), on the other hand, 
equilibria are still spectrally stable, but they are not energy stable. Indeed, the sign of the second variation of the 
free energy functional, in this case, is indefinite because of the presence of NEMs. Equilibria of the second type might 
then be prone to dissipation-induced or nonlinearity-induced instabilities. 

As anticipated in Sec. U one of the advantages of the Hamiltonian formalism for investigating stability and mode 



V. 



CONCLUSIONS 
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signature, is that it is very general and can be applied to any plasma model with a Hamiltonian structure. The method 
described and applied in the present paper, nevertheless, refers only to systems with discrete spectrum. A natural 
and promising future project, would be to analyze mode signature and perform the normal form analysis for more 



40H42I]. Such equilibria 



complex equilibria that support continuous spectrum using the techniques developed in 5 
have more spatial dependence, for example. 
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